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SUMMARY 


An account is given of explicit solutions in terms of Mathieu function 
functions of the problem of two- dimens ional subsonic compressible 
flow past oscillating airfoils. The results are applied to the 
calculation of three-dimensional corrections for the two-dimensional 
theory and the effect of the incorporation of the three-dimensional 
effects on the Mathieu function solution of the two-dimensional 
problem is shown. The developments are formal and must be supple- 
mented by an appreciable amount of numerical calculations before the 
theory can be applied to specific problems. 


INTRODUCTION 


The present report is concerned with the linearized theory of 
oscillating airfoils in two-dimensional and three-dimensional subsonic 
compressible flow. 

The problem of two- dimensional flow was first treated in 1938 by 
Possio (reference 1) who reduced it to an integral equation and obtained 
approx im ate solutions of the integral equation by collocation methods. 
Subsequently, various authors have extended PobsIo’s work and have 
applied other approximate methods to the solution of Possio* s integral 
equation. Detailed references and a survey of the existing results 
may be found in a monograph by Karp, Shu, and Weil (references 2 and 3 ) • 

A second method of treatment deals directly with the boundary-value 
problem of the differential equation for the velocity potential or the 
acceleration potential. If this is done, it is found that an explicit 
solution of the problem may be obtained by introduction of a- suitable 
curvilinear coor dina te system and that this explicit solution Is In 
terms of Mathieu functions. The earliest report containing this 
Mathieu function approach to the problem under consideration appears 
to be that of Sherman and the present author (reference L) . Independently, 
a similar approach, differing in details, was used by Biot (paper pre- 
sented at the Sixth International Congress for Applied Mechanics, Sept. 19^6, 
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"but not available in published form). Tinman (reference 5), and Haskind 
(reference 6). It appears that Haskind has gone furthest in presenting 
the solution in a form suitable for the by no means simple numerical 
evaluation. 

Part of the present report consists in a reproduction, with minor 
modifications, of the existing results in terms of Mathieu functions 
for the problem of two-dimensional flow. One of the reasons for this 
reproduction is the application of these results, in the remainder of 
the present report, to the problem of calculating three-dimensional 
corrections for the results of the two-dimensional theory. 

Just as the problem of two-dimensional flow can be formulated in 
terms of a one- dimensional integral equation, so the problem of three- 
dimensional flow can be formulated as a two-dimensional integral equa- 
tion. This has been done in an earlier report (reference 7). It has 
also been shown in this earlier report that an approximate theory for 
three-dimensional effects can be obtained by making certain approxima- 
tions in the kernel of the integral equation of the three-dimensional 
problem. On the basis of these approximations, the three-dimensional 
problem is reduced to what amounts to a succession of two two-dimensional 
problems. One of these two problems is of the nature of the problem of 
the two-dimensional theory proper. The other problem consists in 
determining the spanwise variation of circulation around the airfoil. 

In the present report It will be shown in which way the Mathieu function 
solution of the two-dimensional problem is affected by the incorporation 
of three-dimensional effects according to reference 7» 

a 

This work was conducted at the Massachusetts Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


X, Y, Z 
H 
t 
U 

u, w 


p o 


Cartesian coordinates 

defined by equation of lifting surface Z = H(X,t) 
time 

main-stream velocity in X-direction 

components of velocity change caused by presence of 
lifting surface 


density of undisturbed fluid 
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V, P 


a 


b 

0 


co 

0 

p 


X, y, z 


M 

k 

* 


1 - m2 
kM 



r = + z2 

*1 

*2 
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pressure and density changes^ respectively., caused by 
presence of lifting surface 

velocity of sound in undisturbed fluid 

semichord 

perturbation velocity potential 
circular frequency of oscillation 
potential amplitude 
pressure amplitude 

dimensionless coordinates defined by equations ( 17 ) 
and (122) 

Mach number of main stream (u/a) 

reduced- frequency parameter (oib/U) 

modified potential amplitude defined by equation (20) 


part of \jr representing noncirculatory portion of 
flow- 

part of -f representing circulatory flow around flat 
plate at rest and at zero angle of attack 

elliptical coordinates 



h 
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g 

Pa 


downvash function defined "by equation ( 4 l) 
pressure amplitude at airfoil 


F, G 
c 

ce m.» ce n 
se mj se n 
■^mn^ ®mn 

^ e m.* Ce n 

^ e m > Se^ 

c cm-> c sm 

®nj ^nu "bn 


functions determined "by equations ( 45 ) 
separation constant 
even, periodic Mathieu functions 
odd, periodic Mathieu functions 

coefficients of Fourier series expressing cej^ and 
se^j, respectively 

1 

Mathieu-Hahkel functions corresponding to and 

ceQ, respectively 

Mathieu-Hankel functions corresponding to se m and 
se n , respectively 

separation constants 

coefficients defined "by equations ( 84 ) and (56) 





mW 

®m 

Sc (l) 

m cm 

J n k)> J n'(n) , 

h, a 
®h 


contribution to pressure at airfoil for nonclrculatory- 
flow component 

noncirculatoiy portion of lift 

coefficient defined by equations (59) a 11 *! (68) 

moment; about midchord 

coefficient defined by equations (62) and (69) 
hinge moment about x = c 
aileron hinge-moment coefficient 

Bessel functions and their derivatives, respectively 
functions defined by equations (70) 
function defined by equation (71) 
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®x 

* 3 ma 


h o <2) - % (2) 


V 2) 

W 

A, B 

am; a n^ Pm; Pn 
Pa^ 2) 

l( 2 ), m(2) 

6 

^0 

X 

Ko, % 
ft 
6 

7* 

n* 

K 


C M 


a, t 


coefficient defined by equation ( 72 ) 
coefficient defined by equation ( 73 ) 
function defined by equation ( 7*0 
coefficient defined by equation ( 75 ) 
arbitrary constant 

Hankel functions of second kind and of zeroth and 
first order, respectively 

coefficients in equation ( 78 ) 

function defined by equation ( 79 ) 

arbitrary parameters defined by equations (1*4-3) 

coefficients defined by equations ( 85 ) 

pressure distribution at airfoil for circulatory- flow 
component 

circulatory portions of lift and moment, respectively 

auxiliary variable of integration 

arbitrary constant in equation (102) 

function defined by equations ( 107 ) and (125) 

modified Bessel functions 

circulation function 

aspect ratio of surface 

dimensionless spanwise coordinate given by equation (123) 

variable of integration 1 

function defined by equation ( 128 ) 

function defined by equation ( 129 ) 

auxiliary variables of integration 
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Q 


G 

P A> P B 
Q A> % 
Qc 


function defined by equation ( 131 ) 
function defined by equation ( 132 ) 
coefficients defined with reference to equation ( 98 ) 
coefficients defined in accordance with equation (ll4) 
coefficient defined by equation ( 139 ) 


R A = 


R B = 


elVCOalce nM 

^Pn e pV cos £ce n (£) 




ll^ 2 ) 


function defined by equation (l46) 


FORMULATION OF TWO-DIMENSIONAL PROBLEM 


Let Z = H(X,t) be the equation of a nearly plane lifting surface 
in tbe path of an inviscid compressible fluid flowing with uniform 
velocity U in the direction of the positive X-axis. Because of the 
presence of the lifting surface, the velocity field (U,0) is changed 
to (U + u,w) . The disturbances caused by the presence of the lifting 
surface are assumed to be small in the sense that the differential 
equations and boundary conditions of the problem may be linearized. 

The differential equations are, in linearized form. 


+ u 

at 

8u _ 

dx ~ 

- -T— ) 

dx \PoJ 

(1) 

^ + u 

8t 

11 

/ol^ 

_ 

dZ \Po/ 

(2) 

+ U §£. 

= -p. 

ysu + 

(3) 

dt 8x 

K ( 

3 \ax dZ/ 


P = £ 

i 2 p 

w 
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In equations (l) to (4), p Q is the density of the undisturbed fluid, 
a is the velocity of sound in the undisturbed fluid, and p and p 
are the changes of pressure and density caused by the presence of the 
lifting surface. 

The boundary condition of no relative normal flow at the lifting 
surface of chord 2b is, in linearized form, 

|X| < b, Z = ±0, v = ^ + U g (5) 

The form of equation ( 5 ) indicates that w may be taken as an even 
function of Z. The differential equations (l) and ( 2 ) then imply that 
u and p may be taken a3 odd functions of Z. As the pressure change 
p must be continuous except when crossing the lifting surface, it 
follows that a further boundary condition may be taken in the form 

• b < |X|, Z = 0, p = 0 (6) 

Further conditions are the conditions of finite trailing-edge 
velocity, 

X = b, Z = 0, u finite (7) 

and the condition that the motion of the lifting surface produces energy 
radiation without reflection at infinity. 

The problem as stated may be solved by means of a perturbation 
velocity potential 0, in terms of ■which 


( 8 ) 


Combination of equations (8), (l), and (2) expresses p in terms of 0 , 
as follows: 



* = - p °(! +u D 


( 9 ) 
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Combination of equations ( 9 ) > (8), (^), and. (3) leads to a differential 
equation for 0 of the form 



Boundary conditions for 0, Besides the radiation condition at Infinity, 
are 



80 - . tt dH 

at = a? ax 


( 11 ) 



^ finite 

8 x 


( 12 ) 


B<|x|, z = o, M+u^ = o (13) 

Attention is now restricted to the case of simple harmonic motion, 
hy setting 

0(X,Z,t) = e iojfc 0(X,Z) (1*0 

Corresponding expressions are written for H and p. 

The differential equation for 0 is 




while the relation Between pressure amplitude p and potential ampli- 
tude 0 is of the form 


p = -p o (ia0 + U g) 


(16) 
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At this stage it is convenient to introduce the following dimen- 
sionless variables and parameters: 






The differential equation (15) then becomes 


(18) 


^ + j 2kM2 80 + j _ Q 

dx 2 dz 2 1 - M 2 8x 1 - M 2 


( 19 ) 


In order to eliminate the first-derivative term in equation (19) a 
function i|r is introduced, defined by 


Setting 




makes equation (19) read 


- 0 

dx 2 Bz^ 


( 20 ) 

( 21 ) 


( 22 ) 


where 


kM 

1 - M 2 


( 23 ) 
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Combination of equations (l 6 ), (20), and (21) gives for the 
pressure amplitude p" the relation 




eipx/iv^f + 

\ OXy 


(24.) 


■where 


v = 


1 - M £ 


(25) 


In terms of the function ijr boundary conditions (ll) to ( 13 ) assune 
the form 


|x| < 1 , z = 


8 z l/l - M 2 ' dx/ 


(26) 


x = 1 , z = . 0 , ^ finite 

ox 


1 < lx I, z = 0 , iVijr + ^ = 0 

dx 


(27) 

(28) 


The condition of no radiation-energy reflection at infinity is 
taken in the form 

z _>+oo^ Tjf ~ f (x, z)e -:i - K:r (29) 

•where r 2 = x 2 + z 2 and -where f tends to zero as z tends to plus or 
minus infinity. 

Determination of the function i|r is facilitated hy introducing 
two functions and such that 

♦ = + * 2 (3°) 

and -where "both ijr^ and ^ satisfy the differential equation ( 22 ). 

The boundary- conditions for \Jr-j_ and ijfg are chosen as follows: 
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x| < 1 


dt-i 

, z = 0, -i = 
dz 

3+2 = 0 

/l - M 2 V ov dz 

(31) 

x = 1, z = 0, 

dt. dtp 

■“• + * finite 

(32) 


dx dx 

1 < |x|, z = 0, 

dtp 

ti = 0 , ivt 2 + — ^ = 0 
dx 

(33) 


It may be seen that the function tp represents the noncirculatory 
portion of the flow to he determined, "while the function tp represents 
the circulatory flow around a flat plate at rest and at zero angle of 
attack. The intensity of the circulatory flow must he chosen such that 
its infinite t railing-edge velocity cancels the infinite trailing-edge 
velocity of the noncirculatory flow. 


BOUNDARY-VALUE PROBLEM WITH REFERENCE TO ELLIPTICAL COORDINATES 


Explicit solutions of the problem as formulated by means of 
equations (22) and ( 30 ) to (33) “ay be obtained through the introduo 
tion of elliptical coordinates I and £. The Cartesian coordi- 
nates (x,z) are related to the coordinates (£,£) by means of the 
equations 


x = cosh £ cos £ 'l 


z = sinh £ sin £ 


J 


(3*) 


The differential equation (22) becomes, in 'terms of the elliptical 
coordinates. 


+ — — + K 2 (cosh 2 £ - cos 2 f)t != 0 
d| 2 d£ 2 


( 35 ) 
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Partial derivatives with respect 
are related as follows: 


to the old and new coordinates 


where 


djr _ dijr t 
5z Sg Sz bt, Sz 

_ dijr S_g <^Tjf b£ 

dx Sg Sx Sg dx 


SI _ cosh I sin g dg _ sinh I cos g 
c)z A 3 Sz A 


SI _ slnh I cos g 
bx A > 


*5 

bx 


-cosh I sin £ 
A 


> 


A = cosh^g - cos^£ 


(36) 


(37) 


The boundary z = 0, | x | ^ 1 becomes the boundary 1=0, 

0 S g = 2 it. The boundary z = 0, x < -1 becomes the boundary 

0 ^ I ■% g = it. The boundary z = 0, 1 < x becomes the boundary 

0^1^“, g = 0. The boundary z = ±®° becomes the boundary I = °°. 
The coordinate curves I = Constant are confocal ellipses with g = 0 
as the limiting ellipse z = 0, |x| < 1. The coordinate curves 

| = Constant are hyperbolas intersecting the ellipses g = Constant 
at right angles. 


Boundary conditions (31) to (33) may now be written in the following 

form 
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The function g(t,) follows from a comparison of equation ( 38 ) and 
equation ( 31 ) in the form 


g(£) = 


Ue" 1 ^ COS S 

Jl - M 2 



m 


From equation (24) it follows that the pressure amplitude p a at 
the airfoil may he written in the form 


Pa = “ 


Po U e ip cos 5 
h 


M (0,£) 


1 *♦«>,» 

sin t, 


C42) 


The problem now is to solve equation ( 35 ) in such a way that 
boundary conditions ( 38 ) to (40) are satisfied. Solutions suitable 
for this purpose will be discussed in the following section. 


MATHIEU FUNCTION SOLUTIONS OF TWO-DIMENSIONAL 
WAVE EQUATION 


In the present section are summarized some known results which will 
be needed in what follows. These, and further developments, not neces- 
sarily employing the same notation as that used here, may be found, for 
instance, in a recent book by McLachlan (reference 8 ). 

Suitable solutions of the differential equation ( 35 ) for ijr are 
obtained by separation of variables, as follows. Setting 


* = F(|)G(£) 


(43) 


there is obtained 


♦ * cosh3| ♦ 2$L 
F(i). G(£) 


K 2 cos 2 £ = 0 


(44) 


Equation (44) implies, with a separation. constant c, the following 
two equations: 

F" + (k 2 cosh 2 | _ c)F = 0 
G" - (k 2 cos 2 £ - c)G = 0 


( 45 ) 
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Possible values of the separation constant c are determined by 
the requirement that the functions G( £ ) are periodic of period 2 it. 

The periodic solutions of the equation deter minin g G may be either 
even or odd functions of *£. They may be written in the following form: 

00 

ce^g) = 21 A mn cos n £ (^ 6 ) 

n=0 


sem(^) = Smn sln 

n=l 


(^ 7 ) 


The function ce m (g) has the same significance for. elliptical 

coor dina tes as has the function cos mg for polar coordinates. The 
function se m (g) is a generalization of sin mg, in the same sense. 

P 

The coefficients A m and Bmn may be written as power series in k . 
The same is true for the separation constants c rTt) and c sm . Evidently 
these statements imply the following limiting behavior of the results: 


ce m (g) = Anna cos mg 
k = 0 ^ se m (g) = Bum sin mg 


m 



The functions ce m and se m may be normalized in various ways. 
The following normalization condition is chosen here: 


jf E-^«U a « - E=m(S)] 2 15 - 1 


(* 9 ) 


The possibility of obtaining an explicit solution of the problem con- 
sidered in this report is due to the following orthogonality properties 
of the Mathieu functions (when m ^ n) ; 



ce n (g) dg = 0 


se n (g) dg = 0 


.(50) 
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C omin g now to the functions F(£), as determined by the first of 
the differential equations (V?)., there will be two linearly independent 
solutions F(£) for each value of the separation constant c. These 
two solutions must be combined in such a way that their behavior at 
infinity is as follows: 


6— F(£) w — e“ iKel ( 51 ) 

\l? ' 

This behavior is required in order that the- energy- radiation condition (29) 
be satisfied by the solution to be obtained. Functions F(|) which have 
this property and which correspond to the functions ce m (£) and Be m (U 
are here designated by Ce m (S) and Se m (|). 

In what follows these two systems of functions are considered 
normalized by the requirement that 

Ce m «(0) = Se m ‘(0) = 1 (52) 

Note that the functions Ce m and Se m may be obtained from func- 
tions Ce^ an d. Se m * which are normalized in a different manner by 
the relations Cem = Ce m yCe m ** (0) and Se m = Se m */se m *' (0). A corre- 
sponding statement applies to the functions ce m and se m . 

There exist various infinite-series representations for Ce m and 
Se m for which reference is made to the literature (reference 8), as no 
use of them is made in the present report. 

As the problem under consideration is linear, all solutions of the 
type here discussed may be combined and written 


t = 5 1 a m ce m(£ )Cem(£ ) + r ^m se m(^)^ e m(^) (53) 

m=0 

The constants a m and bm will be determined by the application of 
appropriate boundary conditions. 


DETERMINATION OF NONCIRCULATORY-FLOW COMPONENT 


Considered first is the simpler problem of determining the func- 
tion which, as can readily be seen, corresponds to a flow without 

circulation. .Boundary conditions (38) to (40) indicate that ijr-^ will 
be an odd function of thus 
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+1 = 21 

m=l 


( 5*0 


Boundary condition ( 38 ) for ijr]_ leads to the following relation: 


21 = s(£) 


( 55 ) 


In view of the orthogonality relations (50) and the normalization 
conditions (49) and (52), the coefficients 'b m are obtained from 
equation (55) in the form 




-r 

J o 


sin £ g(£) se m (£) d£ 


(56) 


The contribution P a ^ to the pressure at the airfoil follows 


from equations (42) and (54) in the following form: 


? ( 1 ) - - £s£ 


^cosfijlv 22 b m se m (£)Se m (0) - 21 >m Be m* (OSeJO) 


(57) 

From equation ( 57 ) it follows that the noncirculatory portion i/ 1 ) 
of the lift may be written as 


l ( 1 > - ab 


/: 

r 


tT ( 1 ) 


dx 


= 2b / P n sin £ d£ 
0 


= -2p 0 U 21 i>m Se m( 0 ) f e iia cos 5 [Iv sin £ se n (£) - se n ‘(^)] 


a£ 


The following abbreviation may be used: 


fc m 


= J e^ cos ^ |iv sin £ se^) - se n *(^r| <i£ 


(58) 


(59) 
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Then, 


L (l) = -^UXVW 0 )^ 


( 60 ) 


A corresponding expression for the moment about midchord becomes 


M^(0) 



— “2p Q Ub b m Se m (0)m m 

where 

% = J cos £ e ip cos sin £ se m (0 - se^ (£)} d£ 
± * l m 


( 61 ) 


( 62 ) 


Analogously, there is obtained for the hinge moment about x = c, 

¥ c ( 1 ) = 2b 2 J (x - c)p a (l) dx = -2p 0 Ub 21 l>mSe m (0)m cm ( 63 ) 

c 

vhere 

r 08 **"^C 

(cos £ - c)e^-^ 008 S|±V sin £ se m (0 - se m *(C)j d£ (64) 

It may be further indicated in which way the coefficients and 

l m may be evaluated. By suitable integration by parts in equation (59) 
there is first obtained 

i m = - j) f e lllc ° s W&> a? (65). 

There are introduced into this equation the series 

se m »(£) = 21 • B mn n cos 
n 


(66) 
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and the relation v/n = M"^. This gives 

I m - (h - Z "W* ** C0S ? C0S n? « (67) 

\M 2 / n J 0 

The integrals in equation (67) are well-known formulas for Bessel func- 
tions and equation (67) may therefore he written in the alternate form 

I* - * - 3) X vAw (®) 

From equation (62) it follows then that 

% = -lj X B m n ni n J n , (^) (69) 

As tables exist for the coefficients and the Bessel functions J n (n) 

and their derivatives J n ’ (n), the coefficients Z m and 1% may 
readily he evaluated for given values of M and k. 

Before a similar statement can he made about the hinge-moment coef- 
ficients mem it will he necessary tc tabulate the function 



e -iH cos £ 


cos m£ d£ 


and its derivative with respect to \i for various values of m and c. 

Having evaluated the coefficients Z m , m^ and m rm , it is next 
necessary to evaluate the remaining coefficients hj^ in equations (60) , 
(6l), and (63), which depend on the motion of the wing, through equa- 
tions (56) and ( 4 l). As will next he shown, this evaluation is readily 
carried out for the two basic degrees of motion 


H(x) = h 
H(x) = abx 



( 70 ) 
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There is obtained 


and 


eh(S) = 


U £^_ 

|/l - m2- 


ikh 


"bmh “ 


Ulkh 
|/l - 


‘•m 


(71) 


(72) 


•where 



e -i^c°Bt sin ? sem (E) at 


“ 51 

n 



e -iii cos £ sin 5 sin n£ d£ 


= - J ^ e -1 ^ cos S [cos (n - l) £ - cos (n + l) Q d£ 






- - § £ £ B mn (-l) n i n J n (^) 


( 73 ) 


Sa(6) “ 


Ue 


-i(JLX 


|/l - M 2 


: ab(ikx + l) 


(74) 


and, with I m from equation (73), 

Uob / 


Slm\ 
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In order to evaluate the corresponding coefficients for aileron 
motion one is again led to nontahulated integrals of the same form as 
those occurring in the evaluation of the aileron hinge-moment 
coefficients m^. 

In summary it may "be said that existing tables are adequate for 
direct calculation of the noncir culatory part of lift and moment for 
bending and torsion motion of the wing section but that additional 
tables are required for aileron motion and for aileron hinge-moment 
calculation. 


DETEEMINATION OF CIECULATOEI'-FLOW COMPONENT 


The essential difficulty in the present problem is the determina- 
tion of the function Boundary conditions ( 38 ) and (40) for i|r 2 

are homogeneous and so the solution ty 2 will contain an arbitrary 
multiplicative constant. This may also be seen by rewriting the boundary 
conditions for ijr 2 in the form 


cH, 


= 0, |x|<l, ^=0 


z = 0, 1 <jx|, ivtyp + = 0 

^ dx 




(76) 


The conditions for 1 < x| may be rewritten in the form 


z = 0, x % -1, t 2 = 0 
z = 0, 1 ^ x, i]r 2 = Ce“^ Vx I 


(77) 


where C is an arbitrary constant and where account has been taken of 
the fact that the uniform main flow is undisturbed far in front of the 
airfoil. 


As ijf 2 is an odd function of z, it follows that ^ and 
cty 2 /c1x are discontinuous across the line (z = 0, 1 ^ x) and this is 
what is responsible for the inconvenience of the problem. 
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It had previously "been proposed (reference 4) to take tyg in the 
following form: 




/ 

J 1 


r ivx' aH 0 (2)(Kr) to . + 2Z^ (2) se n (OSe n (l) 
dz 


( 78 ) 


where 


r 2 = (x - x' ) 2 + z‘ 


The coefficients 1^(2) in equation (78) are then to he determined 

from the remaining boundary condition, which requires that — — = 0 

Sz 

when z = 0, |x| % 1. While this determination is possible in principle, 

it appears to lead to inconvenient formulas when it comes to the evalua- 
tion of the coefficients b^ 2 ). 

An alternate procedure due to Haskind (reference 6) is at present 
believed to be the most convenient approach. Haskind’ s procedure makes 
use of a function W defined by 


Sw 

Sz 


iVtg + 


St 2 

Sx 


( 79 ) 


and required also to satisfy the differential equation V^W + k 2 W = 0. 
As Sw/Sz is an odd function of z, it is indicated that W itself 
is an even function of z. Furthermore, the function W is continuous 
outside the slit z = 0, |x| ^ 1. Consequently, W may be taken in 

the following form: 


W - Y1 V s m( ? ) Ce m(5> 

m 


(80) 


In order to obtain the coefficients ajjj it is observed that the 
first of equations (76) is equivalent to the following boundary condition 
for W: 

z = 0, |x | < 1, = 0 (81) 

' dz 2 
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In view of the form of the differential equation for W, equation ( 8 l) 
is equivalent to the following equation: 

z = 0, |x| < 1, W=A cos kx + B - S ^ Q -- K — ( 82 ) 

K 

Comhination of equations ( 80 ) and ( 82 ) shows that 

amcemd; )Ce m (0) = A cos (k cos 0 + B sin ( 83 ) 

The orthogonality properties of the functions ce m lead then to the 
following expressions for the coefficients a m : 

(84) 


(85) 


W = A X a m ce m ( ^)Ce m ( I) + B 21 Pm ce m( OCeJS) ( 86 ) 


am = Actm + Bp. 


m 


where 


•^m — 


f 

Jo 


cos (k cob £)ce m ((;) d^ 


Cs m (0) 




Pm _ 


[ k *" 1 sin (k cos £)ce m (£) d£ 

J 0 


Ce m (0) 


J 


Now 
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and 


p a ( 2 ) = - ^ eta (±v * 2 + g) 


z=0, |x|<l 


PoU 


^x/dvA 


\dz/ z = 0 , |xj<l 




ip. cos £ | 

f 1 dW\ 


V s in S dV 6=0 

ih cos £ 

k 2Z ct m ce m(^) 

sin £ 


+ B 



Calculation of lift l(^) and moments M^ 2 ^ proceeds in a manner 
which is entirely analogous to the procedure used in equations (58) to 

(75) in order to obtain expressions for and mW. These expres- 

sions may he readily listed as soon as numerical calculations are 
intended. 

The foregoing solution ijrg still contains two arbitrary param- 
eters, A and B, instead of one as it should. The reason for this 
is that, instead of the condition = 0 at the airfoil, only 

the less restrictive condition (iv + d/dx) d^/dz = 0 at the airfoil 

has been satisfied. A relation between A and B which takes account 
of this fact is obtained as follows. From 




5 ^ 5 T " e 

there is obtained by integration 


[Tivx <^2 

L dT 


_x 


= r eta' ax' 

-CO dz 2 


d 2 W 

dz 2 


( 88 ). 


—00 


(89) 
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In view of the form of the differential equation for W, and since 
z)/dz = 0, equation (89) is equivalent to the following equation: 


iVx ^2 _ 
8z 



dx* 


(90) 


Equation (90) , -when suitably applied, leads to the required relation 
between the constants A and B. Since 8^2/8 z = 0 at the airfoil, 

equation (90) implies that 


lim 
z — ^0 u - 


fl eivx ’ + * % ) - °’ M Sl 


(91) 


Since + K% = 0 -when z = 0 and | x I <1, there follows from 

Sx^ 

equation (91) 


lim 
z 0 <J -00 
€ 0 


r 1 +€ e w=!f« + K 

^ -CO \8x 2 


% 


dx = 0 


(92) 


Equation (92) may be transformed by integration by parts in various ways. 
One such form is obtained by writing 



This makes equation (92) appear in the following form: 

(K 2 - v2) f 1 e 1 ” aW(X< 0) to - 
J -ax 

e — iv |a( iVK sin k + k 2 cos k) + B(iV cos ic - k sin k)J 


(93) 


(9^) 
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Another equation which follows from equation ( 9 2) and which is 
equivalent to Haskind’ s equation is obtained by writing 



i Vx 

ax 2 


dx = 


e ivx aw _ iVe iVXy 

ax 


n-i+e 



-00 


e iVx W dx 


(95) 


Combination of equations (95) and (92) gives 
0 


(re 2 - 


V2) J~ e lvy \l(x, 0) 


dx = 


rivll, 


A(iv cos k - k sin jc) - B fc 

L V 


COS K + 


sin 


(96) 


It may be remarked that equation ( 96 ) cannot be used to obtain the 
correct result for the limiting case of incompressible flow for which 
K = 0. The reason for this is, as will be shown, that when k = 0 the 
integral on the left of equation ( 96 ) is not convergent. On the other 
hand, equation (94) does lead to the correct result in the limiting case 
of incompressibility. 


In order to evaluate equations ( 96 ) or (94) it is observed that, 
for £ = «, x = -cosh £ and dx = -sinh 5 d|. Hence, 
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Now take ¥ from equation (86) and obtain* first from equation ( 96 ) * 


A|e“^ v (K sin k - iv cos k) + 

(k 2 - v2) ^T~ aT1 ce n (it) J e" 1V C ° Sh ^Cejj^Cl) sinh I dTj + 

B Je“ iv ^ cos K + sin + 

(k 2 - v2) Pnce n (rt) J e- iVC ° 8 h| Ce n (g) sinh £ dlj = 0 ( 98 ) 

The corresponding result based on equation (9^0 is 

e“^- V K(iv sin k + k cos k) + 

U 2 - V 2 ) <%«*(*) f e- lvc0Bh fee n '(|) di 

d Q 


B 


e -iV(iv cos K _ K sin it) + 


U 2 - V 2 ) 21 Pn“n<*> f e' iV c ° si W (I) « I - 

J o 


iTj = 0 


(99) 


It thus appears that determination of the appropriate values of the 
ratio A/B involves the calculation of the set of infinite integrals 


or 


r e _iv co ® hg Ce n , (l) sinh I a! 
J 0 

r e- iV00BhS Cen , (6) d6 
J 0 


for nl 1 integer values of n and for -whatever values of k and V are 
of interest. 
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It may finally "be remarked that in the work of Haskind . ( ref erence 6) 
it is shown- that equation (98) can he replaced by an alternative equation 
in which, instead of the above infinite integrals, there occur integrals 
between the limits 0, it. The integrand in these alternate integrals 
consists of. products of a Mathieu function and a function which itself 
is defined as an infini te integral involving a Hankel function in the 
integrand. 


DETERMINATION of c irculatory-flow component for 
INCOMPRESSIBLE FLOW, 


The following discussion is intended to show the usefulness of 
Haskind’ s function W in the derivation of the known solution of the 
problem of incompressible flow. 

Setting k = 0, boundary condition (82) becomes 

z = 0, |x| <; 1, W = A + Bx (100) 

or 

|=0, 0 ^ i < 2it, W = A + B cos ( (101) 

Solution (80) becomes, when k = 0, 

w = — (t + io) + f 2 — ®in cos (102) 

* * n=l 

The constant £q arbitrary and may be set equal to 1. Comparison of 
equations (102) and (101) shows that equation (101) is satisfied by setting 


a-L = (n/2)B 
am = 0 


> ( 103 ) 


(m = 2, 3, • 



( 104 ) 


and, consequently, the function W assumes the form 

W = A(l + £) + Be - ^ cos {; 
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(o) 

The pressure distribution p v 1 at the airfoil becomes, according 
to equation (87), 

- (2) = PpU/ 1. 9»\ 

b \sin £ Si4 0 


Po U / A 


(- 


B 


b \ain £ sin £ 


cos £ 


E ) 


Po u / A 




Bx 


x^ /l - x 2 > 


(105) 


Equation (9^-) , 'which gives a relation between A and B, becomes, 
with k = 0, 


Bive-iv = v2 ^ e - lvCosl1 ^£) d£ 


pOO 

= v2 / e" lv COBh ^ (A + Be- 1 ) d| 
J 0 


(106) 


The integrals in equation (106) become, with 


cosh | = X 


at 


= ax/T 


l 2 - 1 


expressible in terms of modified Bessel functions, as follows: 


(107) 


/ 
^ n 


e -iv cosh I d £ 


-/ 


,-ivl 


l/x 2 " 


1 - 1 


dA = Ko(iv) 


(108) 
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dX 


‘ r e ' ±vx (f^r ■ 1 


p-iv 

= Ko(iv) + K^iv) - _ K^(iv) 

P -3V 

= ^(iv) - (IO9) 

Combination of equations (106), (108) , and (109) relates the con- 
stants A and B in the following way: 

BK L (iv) + AK^iv) =0 (110) 

Combination of equations (105) and (110) gives the required explicit 
expression for the pressure distribution at the airfoil for the purely 
circulatory flow component. 

It may finally he remarked that use of equation (96) instead of 
equation (9*0 is not permitted, when K = 0, as the integral on the 
left of equation (96) does hot converge when the function \J of equa- 
tion (10*0 is substituted in it. The reason for this is that W tends 
to infinity for large values of I when k = 0 and does not do so 
when K f 0. 




SATISFACTION OF TEA FLING-EDGE CONDITION 


It remains to obtain one further condition, besides equation (98) 
or equation (99 ) } for the two constants A and B in equation (87) for 

the pressure component Pa^. This remaining condition is the condi- 
tion of finite trailing-edge velocity as expressed by equation (32). 

As i|f]_ and i}^ themselves remain finite at the trailing edge, equation (32) 

may be replaced by the following conditions: 

Ml 

z = 0, x = 1, - — + ivtyp + - — 

ox - ox 


finite 


( 111 ) 
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From equation (5^) it follows that 

1 


z = 0, 


8x 


sin ^ \dt; /|=o 


sin f; 


y~ h n se n »(S)Se n (0) (112) 


From equations (79) and (80) it follows that 

z - 0, iv 2 + sin £ \8 i) | =0 


= — ^ 2Z(Aa n + BPn)ce n (£) (113) 

Comhination of equations (ill), (112), and (113) leads to the following 
further relation between A and B: 

A 2Z c^ce^O) + B 21 P n ce n(°) = JZ V 5 ^ (0)Se n (0) (lib) 

It may he recalled that the coefficients a a and p n are given hy 
equations ( 85 ) and that the coefficients bn are given hy equation ( 56 ), 
where the downwash function g is defined hy equation ( 38 ). 

VALUE OF CIRCULATION JUNCTION 


For the purpose of the subsequent calculation of three-dimensional 
corrections to the two- dim ensional theory, there is needed' the value of 
the circulation function ft, defined hy 


ft = 2e 


iv 



dx = 2e iv i(f(l,0) 


(115) 
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Since ty-]Xl,0) =0 it follows that the circulatory component of flow is 
the only one which contributes to the value of ft and 

ft = 2e iv * 2 (l,0) V ( 116 ) 

From equation (79) it follows that 



(117) 


(118) 


Since ^(-l^O) = 0 * there follows from equation (ll8) the relation 


e 1V * 2 (l,0) = jT e ivcos ^ ( 119 ) 

and, with W from' equation (80), there is obtained for ft 

\ ft = A 21 “n r e 1V COS 5 cen(5) d£ + B ^ Pn f e iV COS ^ce n (C) dt; 

^ J 0 J 0 

( 120 ) 


For numerical applications it will be necessary to evaluate the 
integrals in equation (120) for the chosen values of the parameters V 
and K. This can be done in a manner which is analogous to the method 
of evaluation of the coefficients Z m in equations ( 65 ) to (68). 

CALCULATION OF THREE-DIMENSIONAL CORRECTIONS TO 
TWO-DIMENSIONAL THEORY 


In what follows there is indicated a procedure of incorporating in 
the results of the two-dimensional theory corrections taking account 
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approximately of the three-dimensionality of the flow over wings of finite 
span. This procedure is based on an approximate integral equation of the 
three-dimensional problem as given in an earlier report (reference 7)* 

For the sake of clarity, attention is restricted in -what follows to the 
case of a lifting surface of rectangular plan form, although the results 
of reference 7 were obtained for lifting surfaces with taper and moderate 
sweep as well. 

In considering the three-dimensional problem, besides the coordi- 
nates X and Z for the two-dimensional problem, a spanwise coordinate Y 
is introduced. The plan form of the lifting surface is given by 

|x| $ b ] 

r (121) 

|y| < sb J 

so that s is the value of the aspect ratio of the surface. The following 
two dimensionless coordinates are introduced further: 


and 


y = Jl - m2 1 


( 122 ) 



(123) 


The final result of reference 7 was an integral equation which, for 
the rectangular-plan- form wing, is of the following form: 


g(x,y) 



^(|,y)G(x - 


k) . d| + fl(y)F(x) + 


f 1 

,,-iVx i d2 

J-i an* 


K(y* - q*; k, s, M) dq* 


(12*0 
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The quantities 
G, F } and K 


G(x - 6; k) = 


K(y* - i)*; k, s 


*<>) - jsi r 

z Jq 



St K t a 11 * 3 - V have been defined earlier. The functions X, 
are defined as follows : 


X = 2 


dfr(x,y,0) 

8x 


iitu 

h- 1| 

2 

x - | 




( 125 ) 


(126) 


F(x) = — [ 


e“ iv ^G(x - i; k) d£ 


( 127 ) 


M ) = ly* - n* 

8(1 - m 2 ) y* - n* 


H i 


(2)/ ksM 


tf/l - M 2 


|y* 


- + 


ksM 


— -|y*-n*| 


iks 


|/l - IA C 




ks(y* - t|*) 


[_ /l - M 2 _ 


(128) 


-ia 


| r~ l z I ae~ m \ a2+ Z 2 / 1 \ 

“) d£ + 


zl 


3 -iMl/'i 2 +^ 2 

/r 2 + S 2 


d£ I dT + 


e 


- 1 M\/t 2+ z 2 


t 2 + z 2 


+ iM dT 


T 2 + Z 2 


da 


( 129 ) 
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To effect the solution of the integral equation (124) proceed as 
follows. Write equation (124) in the form 


g 


(x,y) - e _ivx Q = - i l(|,y)G(x - |) 


d| 


(130) 


where 


Q = 



(131) 


G = G - lwe iV F(x) (132) 

Equation (130) indicates that the solution of the three-dimensional 
problem is equivalent to the solution of a two-dimensional problem with 
a modified downwash function g + . Ag, where 


Ag = -e“ ivx Q (133) 

It must next he shown how this modification of the two-dimensional 
downwash affects the previously discussed solution of the two-dimensional 
problem. 

Returning first to the noncirculatory portion of the two-dimensional 
solution, equation (5^), it is evident that the Ag term is responsible 
for a change of the coefficients into + At^. The coefficients bn 

are given as before by equation ( 56 ), except that now they are functions 
of y. 

b n (y) = [ sin £ g(x,y)se n (^) d£ (x = cos £) (134) 

J 0 

The correction coefficients Ab n are correspondingly given by 

Ab n = -Q f sin £ e -iV cos £ sen (£) d£ (135) - 

J 0 

The integrals on the right of equation (135) may "be evaluated once for 
all. The integrals on the right of equation (134) are the same as those 
which must be evaluated for the solution of the two-dimensional problem 
and for which series representations for bending-torsion motions are listed 
in equations (72) and (75)* 
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Coming next to the circulatory portion of the flow, it is evident 
that the only modification here consists in assuming that the arbitrary- 
coefficients A and B in solution ( 86 ) are functions of the span-wise 
coordinate y. The two functions A and B are related to each other, 
as before, in accordance with equation ( 98 ). This relation may be given 
here in the form 

A(y)P A + B(y)P B = 0 (136) 

where P A and P B are defined with reference to equation ( 98 ). 

The next step consists in the satisfaction of the trailing-edge con- 
dition. This condition is analogous to equation (ll4) except that now it 
is necessary to write + Ab n instead of Thus, 

A(y)Q A + B(y)Q B = + Ab n)ae n , (0)Se n (0) (137) 

where Q A and Q B are defined in accordance with equation (llA) . 

In view of equation (135) there may be written, instead of 
equation ( 137 ) , 


A(y)Q A + B(y)Q B + Q(y)Q c = W ( 0 )Se n ( 0 ) ( 138 ) 

where the coefficient Q c is defined by 

Q c = sen* (0)Se n (0) f sin £ e" 1V C ° B 5 se n (£) d£ (139) 

J 0 

The function Q(y) is defined by equation (131). In equation (131) it 
is necessary to express & in terms of A and B, in accordance with 
equation (120). To this end, equation (120) may be written in the form 


Then, 


£2 — AR a + BR b 



dB 


OAo) 


(lAl) 


Q = R A 


K dq* + R b 


-1 dq* 


K dq* 



36 


NACA TN 2363 


Combination of equations (l4l) and (138) gives 


Q A A(y) + Q c R a { K dq* + QgBCy) + Q C R B 



Kdq* = 


bnse n ’ (0)Be n (0) 


(142) 


Equations (142) and (136) represent the main result of the present 
developments, so far as the calculation of three-dimensional corrections 
to the two-dimensional theory is concerned. 

Having calculated A and B in accordance with equations (142) 
and (136), the circulatory component of the pressure distribution follows 
from equation (87), "while the noncirculatory pressure distribution follows 
from equation (57) > "where is replaced by + Ab n . In equation (135) 

for AHq the quantity Q, is to be taken from equation (l4l). 

The foregoing explanation indicates that the major difficulty in 
calculating three-dimensional corrections, once the results of the two- 
dimensional theory are known, consists in evaluating the integrals 


and 



In order to do this, it is necessary to have the numerical values for the 
function H, as defined by equations (128) and (129) • Once the func- 
tion K is calculated, the determination of three-dimensional correc- 
tions to the two-dimensional theory can be effected in a manner analogous 
to what has previously been done for the corresponding problem of incom- 
pressible flow, either in accordance with reference 9 or In accordance 
with any modifications of the scheme used in reference 9» 
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The following further development of the preceding scheme establishes 
a somewhat closer connection with the corresponding results for the. 
problem of incompressible flow. Solve equations ( 136 ) and (l4o) for A 
and B, as follows: 


r a p b - Va 

" P A 

A Q 

r a p b - Va ■ 


(1^3) 


Combination of equations (lV3) and (87) gives the circulatory pressure 
distribution at the airfoil in terms of the circulation function ft. An 
integral equation for the function ft is obtained by combining equa- 
tions ( 138 ) and (1^3 ) , in the following form: 


-k5 ^5^ «(y) + Qo f “ Kin* = h„se„'(0)SenfO) (144) 

Vb’Va J - 1 4,1 


In analogy with the corresponding results for the problem of incompressible 
flow, a function p(k,M) . may be defined by 


n(k,M) = Q c 


R A F B ~ R J^A 

Vb - Va 


(145) 


It may further be written 


fi ( 2 ) (y) ^V'n'C»S%(0) 


(146) 


Equation (l44) then assumes the form 


ft(y) + p(k,M) 



K dq* = ft^ 2 ) (y) 


-1 <h]* 


(1^7) 
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In this form equation (lVf) is of the same appearance as a corresponding 
equation for incompressible flow (reference 9)- Equation (lVf) is valid 
for lifting surfaces with rectangular plan form. The corresponding equa- 
tion for arbitrary plan form, provided sweep is moderate, may also quite 
readily he established on the basis of the results in reference 'J. 


SUMMARY OF RESULTS 


It has been shown that the amplitude of the pressure distribution 
on a rectangular-plan- form airfoil oscillating with frequency oo is 
given by the following formulas: 


Pa 




U) = - ^ e il1 cost IV 21 0>n + Ab n )se I1 ({)Se n (0) - 


sin f; 


+ J Ab n )se n , (^)Se n (0) 


_ (2) Po U cos t 
Pa - “ h 

D sin £ 


3 n ce n^) 


jpji a n ce n^ C ) 

Q(y) 


R A P B " ^A 


In the equations for p ^ and p ^ 


kM 


P = 


M = 


U 


1 - M c 


v = 


1 - M £ 


cos £ = X = — 
b 


v _ a* 
* U 


K = 


kM 


1 - M 
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The Mathieu functions se^ ce n , Se n , and Ce n are defined hy' equa- 
tions (45) to (49 ) , (53 ) i and (54). 

The coefficients afe defined hy 

h n = f g sin £ se n (^) d£ , 

Jo 


where 


and 



Q = 



dn 

*)* 


K dq* 


The coefficients and f3 




are defined hy 

cos (* cos £) ce n '(£) d£ 
Ce n (0) 


P 


n 



sin (k cos 0 ce n (£) d£ 


K Ce n (0) 
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The coefficients P and R in the equation for P a ^ are defined "by 


P A = e“ iv (k sin k - iv cos it) + 


(k 2 - v 2 ) ^ a n ce n (it) f e "^ v cos ^ L ^Ce n (S) sinh I d£ 

^0 


P B = e- ±v 


(cos k + — sin k) 

V K ) 


( k 2 _ v 2) p n ce n (rt) J e“ lv C08lL ^Ce n (|) sinh | d£ 


and 


R a = f e iv C08 ^ce n (ti) d£ 

Rb - 2^ Pn [ e ivcOB ^ce n (0 d£ 
Jn 


1 0 

The function -0 is to he determined from the integral equation 


fl + n(k, M ) f ~~~ K dt]* = O' 


dfl „ , „ _ n ( 2 ) 


J-l 


where 


n(k,M) = Q c 


B A F B ~ ^A 
Q A P B “ % f a 
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in 


and 


n ( 2 ) _ IfXMO ^ V e n '(0)Se n (0) 


and where the coefficients Q^, and Q c are given by 


Q a = y~ a^ce^O) 

Qb = y~ Pn ce n^°^ 

Q c = ^se n « (O)Se n (0) J e“ iV COS £ S e n (g) sin £ d g 


The kernel K of the integral equation is defined by equations (128) 
and (129) and will not be repeated here. 

The results of the two-dimensional theory are contained in the 
foregoing equations. They appear when dfl/dij* = 0 which makes At^ = 0 

in the equation for Pa^ and fl = in the equation for the 

circulation function. 


Massachusetts Institute of Technology 
Cambridge, Mass., May l6, 19^9 
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